Phase diagram of optimal paths 
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We show that choosing appropriate distributions of the randomness, the search for optimal paths 
links diverse problems of disordered media like directed percolation, invasion percolation, directed 
and non-directed spanning polymers. We also introduce a simple and efficient algorithm, which 
solves the d-dimensional model numerically in 0{N^^'''f ^'^) steps where d/ is the fractal dimension 
of the path. Using extensive simulations in two dimensions we identify the phase boundaries of 
the directed polymer universality class. A new strong-disorder phase occurs where the optimum 
paths are self-afBne with parameter-dependent scaling exponents. Furthermore, the phase diagram 
contains directed and non-directed percolation as well as the directed random walk models at specific 
points and lines. 
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The search for optimal paths in a random environment 
is one of the basic problems of statistical physics 1]. The 
task is to find the path between two points in a random 
energy landscape such that the total energy along that 
path is minimal. After a mapping to a graph where the 
edges bear random weights we arrive at a widely studied 
problem of discrete mathematics (see, e.g., From the 
point of view of physics, the interest comes from intimate 
relations to various fields including thegeometry of do- 
main walls in disordered magnets y, 0, U i vortices in su- 
perconductors y or rupture lines |3- For d-dimensional 
directed lines without overhangs, often referred to as "di- 
rected polymers in random media" (DPRM), there is an 
exact mapping to stochastic surface growth [Sj, enabling 
further connections to fractal surfaces |^, to driven par- 
ticle s_TOtems ^3 to the stochastic Burgers equa- 
tion The main questions for all these applications 
are: How can the geometry of the optimal paths be de- 
scribed and how do the energy fluctuations scale with 
their length? Some realizations of optimal paths for dif- 
ferent disorders are shown in Fig. ^ 

The situation is quite clear for the DPRMs: The 
line is a self-affine fractal resulting in scaling endpoint- 
fluctuations. Imagine that the polymer is held fixed at 
one end and we ask for conformations of length L over 
different disorder realizations. Let (Acc)^ be the aver- 
age squared endpoint-fluctuations of the polymer, then 
Ax oc where Q is the roughness exponent. Similarly, 
for the average squared energy fluctuations (A£')^ we 
have AiJ cx L'^. In d = 2, the exponents are known 
exactly (dprm = 2/3 and uj = 1/3. In higher di- 
mensions numerically calculated values are available |12| . 

Recently, the question of the conformation of span- 
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FIG. 1: Optimal paths for different values of the disorder 
parameters S and e, defined in Eq. (3). 



ning non-directed polymers has attracted considerable 
interest. Indeed, the original optimal path problem does 
not exclude any geometry a priori, in principle overhangs 
should be allowed. In fact, it can be easily shown that 
overhangs have to occur on a sufficiently long line. How- 
ever, recent numerical work on such non-directed poly- 
mers in a random medium (NDPRMs) have given evi- 
dence that, for bounded unimodal distributions of dis- 
order, these overhangs are irrelevant |0 and the ND- 
PRMs are in the universality class of DPRMs [H 
The latter conclusion was also drawn from a real space 
renormalization group study (13.] . 

Most of the studies of optimum paths have been con- 
fined to Gaussian distribution of the randomness. The 
effect of changing the distribution has been studied in 
some recent papers. In spite of earlier reports 21], uni- 
versality of the DPRMs was found for wide family of 
bounded distributions It was shown JJJ that, when 
the distribution has a power law decaying part for neg- 
ative energies characterized by an exponent fi, the effect 
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FIG. 2: Phase diagram in the disorder parameter S and the 
directedness parameter e, defined in Eq. (3). DPRM: Directed 
Polymer in a Random Medium, IP: Invasion Percolation, DP: 
Directed Percolation, NUOP: Non-Universal Optimal Paths. 



described by Zhang [20j for surface growth can be ob- 
served for DPRMs too: The scaling exponents of the 
paths depend continuously on fi in an intervall of /i. For 
a bimodal distribution, where the energy can be either 1 
or (with probability p and 1—p, respectively) there are 
crossov er p henomena near to the percolation threshold 
Pc |l3l |23| |. Long range correlations in the randomness 
may lead to nonuniversal behavior . 

The standard models of disorder, directed percolation 
(DP) or ordinary percolation can also be formulated 
in terms of global optimization problems , as long as 
the spannin g pa ths are concerned. We briefly summarize 
these cases The optimal path V in the polymer 

problem is given by 



mm 

V 



(1) 



while in the percolation problems we have for the span- 
ning paths P 



min max 1 Ei , 

V \ iev J 



(2) 



where for the directed versions the condition of directed- 
ness to (1) and (2) have to be superimposed. The only 
difference between these equations is that in (2) we have 
the maximum instead of the sum in (1). When looking 
for the bottlenecks along each portion of the path we 
obtain exactly an optimal percolation path. However, 
this path does not scale like the shortest (or chemical) 
distance on a critical percolation cluster, but rather as 
one on an invasion percolation (IP) tree [sl. l26j. or trap- 
ping invasion percolation |26j | . These invasion percolation 
problems have much in common with ordinary percola- 
tion, e.g., the critical point and the fractal dimension of 
the critical cluster agree. 

The geometry of the optimal percolation paths is also 
known. For DP we have again a self-affine line with 
Cdp = 0.633 in two dimensions calculated from the 
ratio of the perpendicular and the parallel correlation 



length exponents. For IP we have a fractal line with a di- 
mension df > 1. In two dimensions its value is df = 1.22 
Ull^. Naturally, the roughness exponent for a fractal 
line is Qp = 1. 

It is tempting to try to formulate an unified picture 
of these diverse but obviously related phenomena and in 
this Letter we address this point. Furthermore we would 
like to explore the possibility of new universality classes 
and at the same time to check earlier results by improving 
the numerics. 

Let us look for a parametrized distribution of random- 
ness such that both directedness and strong disorder nat- 
urally occur as limits. In order to achieve this, weights 
are assigned to each bond as 

E = r^ +e , (3) 

where r G (0,1) is taken from an uniform random dis- 
tribution; £ > and — 00 < (5 < 00 are the directedness 
and the disorder parameters, respectively. This choice 
corresponds to a PDF in the reduced variable E = E — e: 



p{E) 



(4) 



Here E e (0, 1] for 6 > Q while E > I lov 5 < 0. The 
optimal path is searched for in a random environment 
characterized by the bond weights in (3), resulting in 
complex phase diagram, Fig. |21 

Already a simple analysis of (3) reveals interesting con- 
nections to different models of statistical physics. Let us 
take first 5 — 0. Along this line a penalty is to be paid for 
any overhang. However, all directed paths have the same 
weights. This is exactly the statistics of directed random 
walks (DRW) , leading to self affine lines with roughness 
exponent (^drw = 1/2 [2^ . 

The case of NDPRM with a bounded distribution can 
be identified with e = 0, < (5 < 00. The optimal 
paths are in the universality class of DPRM [H 
Obviously, this is even more so when e, the parameter 
suppressing overhangs, is switched in. Interestingly, no 
crossover phenomena at 5 — > are to be expected since 
however small 5 is the global optimization over the ran- 
dom weights matters, since the e weights represent only 
a constant background. 

The limit (5 — > +00 the distribution is more and more 
shifted towards a delta function at 0. Nevertheless, this 
limit corresponds to strong disorder as it can be shown 
with reference to the order statistics |2^ : The cumulative 
distribution is P{E) = p{E')dE' = E^l^ (we first 
assume e = 0). We generate an ensemble of groups of 
-E-values, N in each group. We order the numbers in 
each group in an ascending order so that < i?(2) < 
• • • < -E(jv) • The average value of the nth element in 
this ordering is (£'(„)) = [n/(A^+ 1)]''. Hence, the ratio 
between consequtive elements in this ordering is given by 
{Ef^n+i)) I {E(n)) = (1 + ^/nY which diverges for any n 
as 5 —>■ +00. Consequently, any sum of elements drawn 
from this distribution will be dominated by the largest 
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FIG. 3: Ax as a function of system size L. The slopes of the 
(5 = —2, e = 0) and (—3, 0) hnes are 0.67, while the slopes of 
the {5 = —5, e = 0) and (—6,0) lines are 0.75. 



element, showning that we are in the infinite disorder 
limit, and the optimum paths are percolation paths [l6l| . 
For e = it is the invasion percolation shortest path 
|26j | . However small value of e is taken it will immediately 
suppress the overhangs transforming the optimal path to 
a critical directed percolation path. 

In the regime 5 < the distribution has a power law 
tail with the consequence that only the k < l/|(5|-th 
moments remain finite. If this property is inherited to 
the optimal paths, this would mean that there should be 
nonuniversal behavior as a function of S since the univer- 
sality classes are characterized not only by the exponents 
but also by the universal scaling functions. Moreover, the 
large number of huge obstacles generated by the tail of 
the distribution may lead to an enhancement of over- 
hangs resulting in new geomtries. We apply numerical 
tools in order to see clearly in this point. 

Before that, let us summarize the behavior ai 6 — > —oo. 
In this limit, the distribution becomes not normalizable 
indicating the dominance of large values of E, i.e., this 
corresponds again to the infinite disorder limit, where 
the arguments of order statistics apply straitforwardly 
and we arrive at (2) which is the percolation case. It is 
a question of interest whether the paramter e is relevant 
for this limit and we will come later back to this point. 

Let us turn to our two-dimensional numerical simula- 
tions. The obvious difficulty in the determination of op- 
timal paths is in the global nature of the problem. Given 
an optimal path of span L, the solution for size L + 1 
can be entirely different. However, there are algorithms 
which calculate the optimal path on graph with random 
positiv weights on the edges in polynomial time Qj- We 
have developed an efficient algorithm which is adequate 
to the geometry we considered. 

We now describe this algorithm in some detail. It is 
closely related to the algorithm described by Hansen and 
Hinrichsen in for determining the non-directed per- 
colation threshold. We describe the algorithm for the 
two-dimensional case, i.e., for a 45° tilted square lattice. 
The sample has a cylindrical geometry (periodic bound- 



FIG. 4: AE as a function of system size L. The slopes are for 
(5 = 2,£ = 0) 0.33, for (-1/2,0), 0.25 and for (-1,0) 0.13. 

ary conditions in the horizontal direction). Each bond 
i is assigned an energy Ei. We associate a variable Va 
with each node a. Initially, these variables are all set to 
zero with the exception of the nodes forming the hori- 
zontal edges of the lattice. These nodes are assigned a 
value which is larger than any energy Ei appearing in 
the system. Furthermore, the node forming the anchor- 
point of the polymer, a = is assigned a negative value 
Vb = vo = —\vo\. This node sits on the lower edge of 
the lattice. We then iteratevely update the internal ver- 
tices (i.e. those not being on the horizontal edges of the 
lattice) according to the scheme 

Va min(Va(i) -|- i?a(l), ^q(2) + Ea{2), 
Ki(3) + Ea(3),Va(4:) + i?Q(4)) , 

where a(l) to a(4) are the four nodes that neighbor node 
a and the four bonds joining them. The updating pro- 
ceeds until the node variables no longer change. The 
value of each node variable is then equal to the energy of 
the optimal path between that particular node and the 
anchor node plus the value vq. The number of iterations 
goes as Lx the length of the shortest path, i.e., as L^+'^-f , 
where d/ is the fractal dimension of that path. 

In Fig. 3 we show end point fluctuations Ax as a func- 
tion of system size L for a number of different disorders. 
For 6 > —2 and e = 0, we find Ax ~ L"-^^, indicat- 
ing self affinity with the same roughness exponent as in 
the DPRM problem, where C,dprm = 2/3. We also find 
the same behavior with 5 = —8 and e — 10"*. How- 
ever, for disorders with J < — 5 or less and £ = 0, we 
find Ax ^ L" ''^, which indicates self afflne behavior, 
but with a new roughness exponent which is significantly 
larger than the DPRM one. The self affinity of the op- 
timal paths is further supported by monitoring the av- 
erage length of the paths, I as a function of system size 
L. For all finite disorders investigated, we find I ~ L^ ™. 
Hence, the curve is not a self-similar fractal, and the 
overhangs are irrelevant. Fig. 1 shows different optimal 
paths with different disorders. Even in the most extreme 
cases {6 = —8 and e = 0), there are few overhangs. 

The fluctuations in energy, AE, scales as AE ^ L^^^ 
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in the DPRM case. For 5 > -1 we find that AE is 
self averaging. In Fig. 4, we show AE as a function of 
L for (5 = 2, (5 = -1/2, and 6 = -1, aU with e = 0. 
For S — 2 case, we find an exponent consistent with the 
DPRM value, 1/3. However, for the 5 = -1/2, we find 
a value which is much smaller, 0.25, while for (5 = 1, we 
find 0.13. This may be a crossover effect. 

For values of i5 < —1, AE is not self averaging, in- 
dicating no universal scaling functions exist. However, 
the scaling of Aa:; with two different values of scaling ex- 
ponent suggests that there are still partial universality 
in the system for small values of S and e = 0, with a 
roughness exponent of 0.75. 

Universality classes are characterized by scaling func- 
tions, not only by exponents '2^. Therefore we extended 
our simulations to higher moments of the endpoint fluctu- 
ations (wherever it was possible), because the amplitude 
ratios are known with high precision ^231 and also to the 
study of the moments of the energy fluctuations. From 
our studies we conclude: i) For < S < oo the amplitude 
ratios are within the numerical accuracy the same as for 



DPRM, conflrming universality, ii) For > 6 > — oo 
the fc-th moments of the energy fluctuations become di- 
vergent for k > l/\S\ (i.e., even the second moment is 
nonexistent for 6 < 2) indicating nonuniversal behavior 
(see. Fig. 121 Non-universal optimal paths, NUOP). The 
measurements were carried out for e = 0. In the case of 
e > the situation is even worse: While e does not con- 
tribute to the energy fluctuations, its effect is enhanced 
directedness, i.e., the path is forced through higher E 
valued regions. However, it seems that unversality in 
a restricted sense is present: The roughness exponent 
seems to agree with (dprm for (5 > — 4 and £ = 0. Also 
for large values of e the C becomes close to 2/3 even for 
large negative values of S. Furthernmore, the roughness 
exponent for 6 < —4 is ~ 0.75. Further investigations 
are needed to explore fully this restricted universality. 
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